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Abstract. We have analyzed the far-ultraviolet spectrum of two Galactic 04 stars, the 04If+ supergiant HD 190429 A and 
the 04V((f)) dwarf HD 96715, using archival FUSE and IUE data. We have conducted a quantitative analysis using the two 
NLTE model atmosphere and wind codes, Tlusty and CMFGEN, which incorporate a detailed treatment of NLTE metal line 
blanketing. From the far-UV spectrum, we have derived the stellar and wind parameters and the surface composition of the two 
stars. The surface of HD 190429A has a composition typical of an evolved O supergiant (nitrogen-rich, carbon and oxygen- 
poor), while HD 96715 exhibits surface nitrogen enhancement similar to the enrichment found in SMC O dwarfs which has 
been attributed to rotationally-induced mixing. Following studies of Magellanic Cloud O stars, we find that homogeneous wind 
models could not match the observed profile of O V/11371 and require very low phosphorus abundance to fit the P V/LU118- 
1 128 resonance lines. We show, on the other hand, that we are able to match the O v and P v lines using clumped wind models. 
In addition to these lines, we find that N iV/11718 is also sensitive to wind clumping. For both stars, we have calculated clumped 
wind models that match well all these lines from different species and that remain consistent with Ho- data. In particular, we 
have achieved an excellent match of the P v resonance doublet, indicating that our physical description of clumping is adequate. 
These fits therefore provide a coherent and thus much stronger evidence of wind clumping in O stars than earlier claims. We 
show that the success of the clumped wind models in matching these lines results from increased recombination in the clumps, 
hence from a better description of the wind ionization structure. We find that the wind of these two stars is highly clumped, as 
expressed by very small volume filling factors, namely /„ = 0.04 for HD 190429A and /«, = 0.02 for HD 96715. In agreement 
with our analysis of SMC stars, clumping starts deep in the wind, just above the sonic point. The most crucial consequence of 
our analysis is that the mass loss rates of O stars need to be revised downward significantly, by a factor of 3 and more. These 
lower mass loss rates will affect substantially the evolution of massive stars. Accounting for wind clumping is essential when 
determining the wind properties of O stars. Our study therefore calls for a fundamental revision in our understanding of mass 
loss and of O-type star stellar winds. 

Key words. Stars: winds - Stars: atmospheres - Stars: early-type - Stars: fundamental parameters - Stars: individual: 
HD 190429 A, HD 96715 



1. Introduction 

The spectrum of hot, massive stars reveals signatures of dense, 
highly-supersonic mass outflows that are driven by the strong 
radiation field. These stellar winds have a significant effect on 
the evolution of massive stars, as O stars will lose a sizable 
fraction of their total mass during their lifetime. Mass loss rates 
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are, therefore, a crucial parameter of stellar evolution models. 
These mass loss rates, like other basic stellar parameters, can 
be derived from spectrophotometric analyses based on model 
stellar atmosp heres. For a recent review of the properties of hot 
star winds, see iKudritzki & Pulsl fcOOCh . 

During the last decade, major advances have been ac- 
complished in constructing more realistic model atmo- 
spheres, which consistently account for the stellar wind, 
for departures from the Local Thermodynamic Equilibrium 
(LTE) , and for petal line blanketing ( Hubenv & Land 1 1995: 
[ Hillier & Millerlll998t iHauschildt et alJll997i IPauldrach et alJ 
2001; Koesterke et al. 2002; Hube nv et alT l2003t) . Current 
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state-of-the-art non-LTE (NLTE) model atmospheres assume 
smooth, homogeneous and stationary winds. There are, how- 
ever, theoretical arguments as well as observational evidences 
that the winds of hot stars are extensively structured. Radiation 
hydrodynamics simulations of the line-driven flow's nonlinear 
evolution show the line force to be highly unstable and lead 
to strong reverse shocks ("see. e.g.. | Ow ocki et alll988l lOwockil 
1994lFeldmeierlll995UOwocki &Pulslll999l) . We thus expect 
the presence of density contrasts in the wind of O stars, which 
is supported by several lines of observations. 

First, the soft X-ray emission of O stars is widely believed 
to originate f rom shocks propagating t hrough the stellar wind 
jLucvlll982l ICassinelli & Swank Il983l) . With the advent of 
high resolution X-ray spectroscopic capabilities on Chandra 
and XMM-Newton, detailed wind-shock models can be tested 
by fitting X-ray line profiles. The fitted lines indicate that the 
hot emitting plasma is lo cated th roughout the wind starting 
close to the photosphere dKramer et alJl20 03). However, the 
high opacity of their wind mo del has to be lowered in order to 
reproduce the X-ray flux level. lOskinova et all J2004) recently 
offered a solution to this issue considering highly fragmented 
winds. 

Second, temporal variability of O-star wind line profiles 

is well documented and might indicate the propagation of 

disturbances throughout the wind, most ofte n observed as 

Discrete Absorption Co mponents (see, e.g., IHowarth et alJ 
. , , , , _i i . . , 1 1 . , , , , _i 



1995; Massa et al. 1995). In the spectrum of the prototypical 
O supergiant £ Pup. lEversberg et al] |l998) observed stochas- 
tic variable substructures in He n/14686 that move away from 
the line center with time. They explained them as evidence of 
blobs or cl umps moving outward in the stellar wind. 

Third, Sillier (1991) showed that the weak electron scat- 
tering wings of He n/14686 in the spectrum of Wolf-Rayet 
(WR) stars can be explained by clumped wind models and 
can thus be used to diagnose density inhomogeneities in hot 
star winds. The inclusion of clumping in wind studies of WR 
stars resu lted in lowering their mass loss rates by a factor 
of 2 to 4 jMoffat & Rober1ll994llHamann & Koesterkdl 19981 
iHillier & MilleJl999l) . 

Lastly, spectral evidence of clumping has been re- 
ported in recen t studies of Magellan ic C l oud (MC) O stars. 
ICrowther et all (H5o2b) and IHillier et alJ J2003h found that 
phosphorus needs to be significantly underabundant relative 
to other heavy elements in order to match Pv/L11118, 1128 
in O supergiants. Alternatively, theoretical line profiles from 
clumped wind models match the observed P v lines, assum- 
ing a "standard" phosphorus abundance. ^From studyin g the 
wind i onization of a large sample of LMC O stars. iMassa et alJ 



(2003) came to a similar conclusion about the P v resonance 
lines provid i ng an indication of wind clumping. Moreover, 
IHillier etalJ J2003) showed that the clumped wind mod- 
els match strong Fe lines, while homogeneous wind mod- 
els predict marke d blue asymmetr ies that are not observed. 
Independently, iBouret et alJ J2003I) successfully matched for 
the first time the O V/11371 line profile using clumped models. 
A good match to this line was never achieved with homoge- 
neous wind models. Bouret et al. found that clumping starts just 
above the photosphere. They derived mass loss rates for several 



O dwarfs in the SMC cluster NGC 346, which are a factor of 3 
to 10 smaller than those obtained from smooth winds. This re- 
duced mass loss is expected to significantly alter the predicted 
ev olution of these sta rs . 

|PulsetalJ in996) showed the existence of a tight correla- 
tion between the modified wind momentum and the stellar lu- 
minosity. This Wind momentum - Luminosity Relation (WLR) 
varies for Galactic and MC O stars because of the different 
metallicities. Stars with lower abundanc es of heavy elements 
have weaker winds ( Walbor n et alJl200 0). Mass loss rates pre- 
dicted using radiative line-d riven wind theory bring support to 
the WLR ( Vin k et alj|20 00). hence allowing a derivation of O 
star luminosities and of extragalactic distances. However, the 
dependence of the WLR with the luminosity class remains an 
open problem. Recent theoretical studies do not predict such a 
dependence while, empirically, a differe nce in the WLR for su- 
pergiants and for dwarfs has been found. Repol ust et alj j2004) 
argued that this discrepancy between the empirical and theo- 
retical WLR's might be solved if the mass loss rates derived 
from Ha are affected by clumping in the lower wind region. 
Clumping would result in decreasing the mass loss rates of 
Galactic O stars with Ha emission by a factor of about V5. 

Despite mounting evidence that hot star winds do not ap- 
pear to be smooth and homogeneous, most studies still assume 
so with the hope of determining average wind properties at 
least. In MC O stars, however, spectral analyses accounting for 
clumping yield significantly lower mass loss rates, with poten- 
tially crucial consequences for our understanding of massive 
star evolution, and for the radiative line-driven wind theory and 
its application to extragalactic studies using the WLR. In this 
context, we believe that it is of special importance to extend 
the recent studies of MC O stars to Galactic O stars which have 
stronger winds. Furthermore, it is essential to show that a con- 
sistent picture may be drawn from the different diagnoses of 
clumping, in particular P V/l/11118, 1128, O V/11371, and Ha, 
ensuring that the spectral diagnoses of wind clumping are not 
spuriously affected by abundance or ionization effects. 

In this paper, we investigate the properties of two Galactic 
04 stars, HD 190429A and HD 96715, in order to derive their 
properties and address these issues. These two stars were cho- 
sen to cover a wide range of luminosity classes, thus investi- 
gating wind clumping in O dwarfs and in O supergiants, and to 
compare Galactic stars to SMC 04V stars for which Bouret et 
al. found evidence of clumping. The general stellar properties, 
along with the observational data and reduction are presented in 
13 Sect. [5] considers our model atmospheres. We discuss the 
methodology used to derive the stellar parameters, chemical 
abundances, and wind parameters, and the related uncertain- 
ties, for each star, in ^ Our results are then put into a broader 
context in |5] comparing in particular the derived mass loss 
rates to other predictions based on homogeneous wind models. 
General conclusions are summarized in ^6l 



2. Stellar sample and observations 

We have selected two Galactic 04 stars, HD 190429 A and 
HD 96715, which have been previously observed in the far 
ultraviolet (FUV) by FUSE and WE, to investigate the FUV 
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Star 


RA 


Dec 


Spectral 


FUSE 


V 


B-V 


E(B - V) 


logA^H i) 


logW(H 2 ) 




(J2000) 


(J2000) 


Type 


ID 


(mag) 


(mag) 


(mag) 


(cm" 2 ) 


(cm -2 ) 


HD 96715 


11 07 32.9 


-59 57 49 


04V((f)) 


PI 024301 


8.27 


0.10 


0.42 


21.3 


19.7 


HD 190429 A 


20 03 29.4 


+36 01 30 


04If+ 


PI 028401 


7.12 


0.14 


0.51 


21.5 


20.3 



spectral signatures of clumping (O v, P v). The two stars are 
hot enough (04 spectral type) to reveal O V/11371 in their spec- 
trum. We have selected one supergiant and one dwarf pro- 
viding an initial comparison between the two classes. The 
two stars are part of the FUSE atlas of Galactic O stars 
JPeilerinetal.1120021) . from which we extracted basic stellar 
data (Table[D- HD 190429A is a bright O supergiant for which 
extensive o bservational and modeling work has been pub- 
lished (e.g., IWalborn & HowarfhlbOOOt iMarkova et al.ll2004l 
iGarcia & Bianchil2004 . lGarcia & Bianchi d2004 recently re- 
ported the first detailed analysis of the 04V((f)) star HD 96715. 

We have extracted IUE short wavelength, high resolution 
spectra from the Multimission Archive at the Space Telescope 
Science Institute (MAST). The SWP spectra cover the spec- 
tral range, /Ml 150-2000 A, at a resolving power R = 10,000. 
Each spectrum was processed with the NEWSIPS package. We 
have selected all spectra obtained through the large aperture, 
that is, 16 spectra for HD 190429A and 4 spectra for HD 96715 
(Table[2j. These spectra do not show conspicuous signs of vari- 
ability, in flux as well as in line profiles. Although we can- 
not rule out definitively variability (in particular for the super- 
giant HD 190429 A), this absence of variations justifies our co- 
adding all merged extracted spectra so as to form an average 
spectrum for each star. Data points flagged by the NEWSIPS 
software have been excluded; in particular, this concerns the 
saturated portion of the long-exposed images SWP 43980 and 
SWP 43981 at A Z 1700 A. Finally, we smoothed the co-added 
spectra to a resolution of 40 kms -1 in order to increase the 
signal-to-noise ratio. The spectra of both stars show a large 
number of narrow lines of interstellar (IS) origin. This IS con- 
tamination becomes an issue especially for HD 96715 because 
its apparent rotational velocity is relatively low (vsin i= 80 
kms~ 'i lHowarfh & Priniall989l) . 

The processed FUSE spectra have been retrieved from 
MAST too. The nominal spectral resolution is 20,000, or about 
20 kms -1 . Details about the observations and data reduc- 
tion can be f ound in the FUSE a tlas of Galactic OB spectra 
presented by Pell erin et alJ (120021) . Individual sub-exposures 
have been co-added for each segment and then merged to 
form a single spectrum, using Lindler's FUSE-REGISTER pro- 
gra m. We avo ided contamination by the so-called worm arte- 
fact JSahnow et al. 2000) by using only the LiF2A spectra on 
the long-wavelength side (/U1086-1 183 A) of the spectrum. 
Finally, the co-added merged spectra have been smoothed to 
a 30 km s _1 resolution to enhance the signal-to-noise ratio. 

^From the IUE and FUSE spectra, we have measured the 
atomic and molecular hydrogen column densities towards the 
two stars, fitting Lya and H2 lines, respectively. The values are 
listed in the rightmost columns of Table ^ indicating that the 
FUSE spectra suffer from severe blending from the broad ab- 



sorption bands of H2. This absorption especially becomes a se- 
rious issue shortwards of 1006 A (SiC2A channel), and thus 
little information about the stars can be inferred from the anal- 
ysis of this part of the spectra. Furthermore, numerous IS metal 
lines are also present in the FUSE spectra, superimposed to 
the stellar lines. These lines are narrower than t he stellar lines 
but, b ecause of their large number (see Table 3, Pel lerin et alJ 
120021) . they may complicate the analysis of the stellar compo- 
nents (especially for weak lines in HD 96715). 

3. Modeling assumptions 

Unified models, with a consistent treatment of the photosphere 
and the wind, are mandatory to analyze the P Cygni profile of 
strong lines in the UV spectra of O stars and, thus, to derive 
the basic wind parameters (mass loss rate, terminal velocity). 
However, the bulk of spectral lines in O stars are formed in the 
photosphere where velocities are small and geometrical exten- 
sion is negligible. Hydrostatic, fully-blanketed, NLTE photo- 
spheric models may then remain a preferable alternative to uni- 
fied models with a simplified treatment of metal line blanket- 
ing. We demonstrated this point in the case of low metallicity 
stars in the SMC, where stellar winds are weaker than those of 
Galactic stars JBouret et all2003 j. Here, we deal with Galactic 
objects and, in particular, with an extreme 04If+ supergiant 
(HD 190429 A) which exhibits a strong stellar wind. This is 

Table 2. IUE SWP spectra. 



Star 


Image 


fexp 


Obs. Date 






(s) 




HD 96715 


SWP 21999 


3000 


1984-01-13 




SWP 22000 


1800 


1984-01-13 




SWP 43980 


9000 


1992-02-13 




SWP 43981 


9000 


1992-02-13 


HD 190429 A 


SWP 04903 


1515 


1979-04-09 




SWP 38958 


1500 


1990-06-02 




SWP 38965 


1500 


1990-06-02 




SWP 38970 


1500 


1990-06-03 




SWP 38973 


1500 


1990-06-03 




SWP 38978 


1500 


1990-06-03 




SWP 38981 


1500 


1990-06-03 




SWP 38986 


1500 


1990-06-03 




SWP 38989 


1500 


1990-06-03 




SWP 38994 


1500 


1990-06-04 




SWP 38998 


1500 


1990-06-04 




SWP 39003 


1500 


1990-06-04 




SWP 39006 


1500 


1990-06-04 




SWP 39011 


1500 


1990-06-04 




SWP 54573 


1500 


1995-05-02 




SWP 55986 


1500 


1995-09-22 
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thus an exemplar case for studying the wind contribution to 
weak photospheric lines, comparing the predictions of photo- 
spheric and unified models to observations. We have performed 
such an analysis using model atmospheres calculate d with 
the photospheric program, Tlusty (Hubenv&Lanz 1995), 
and w ith the unified model code, CMFGEN farther & M iller 
1998). The atomic data used by the two codes mostly come 
from the same sources, thus ensuring a meaningful comparison 
of the stellar parameters derived by the two codes. CMFGEN 
does not solve the full hydrodynamics, but rather assumes the 
density structure. We use a hydrostatic density structure com- 
puted with Tlusty in the deep layers, and the wind part is de- 
scribed with a standard y6-velocity law. The two parts are con- 
nected below the sonic point at v(r) » lSkms" 1 . For mor e 
detai ls on these two codes, we refer to Hubenv & Lanz ( 1995), 
Lanz & Hubenvl J2003I) IHillier & Milled dl998h . Iffillier et alJ 
J2003I) . and lBouret et alJfe003l) . 

Radiatively driven winds are intrinsically subject to insta- 
bilities, resulting in the formation of the discrete structures 
or "clumps". As discussed in §1, there are several observa- 
tional evidences as well as theoretical arguments that foster 
the concept of highly-structured winds. To investigate spec- 
tral signatures of clumping in the winds of HD 96715 and 
HD 190429 A, and its insuing consequences on the derived 
wind parameters, we have thus constructed clumped wind mod- 
els with CMFGEN. A simple, parametric treatment of wind 
clumping is implemented in CMFGEN, which is expressed by 
a volume filling factor, /, and which assumes a void interclump 
medium and the clumps to be small compared to the photons 
mean free path. The filling factor is such that p = fp, where p 
is the homogeneous (unclumped) wind density. The filling fac- 
tor decreases exponentially with increasing radius (or, equiva- 
lently, with increasing velocity) 

/ = foo + (1 - /ex.) exp(-v/Vd), (1) 

where v c i i s the velocity at wh ich clumping starts. IHillier et alJ 
( 2003| and lBouret et at) |2003 ) found that clumping starts close 
to the photosphere. We have thus adopted, v c i = 30km s _1 , that 
is, clumps start forming just above the sonic point. 

The observation of strong resonance lines of highly-ionized 
species, like O vlM1032, 1038 and N V/U1238, 1242, provides 
an indirect evidence of X-ray emission in the stellar winds. 
An important consequence of the X-ray and EUV shock radi- 
ation is to enhance photoionization that results in "wind su- 
perionization". Auger proc esses are accounted for i n calcu- 
lating the wind ionization llCassinelli & O lson 1979). X-rays 
probably originate from the cooling zone of post-shock re- 
gions, while shocks aris e from radiative instabilities inher- 
ent to line-driven w inds JLucv & Solomonlll970t IlucvII1982|: 
IOwockietal¥ l988). X-ray emission has been observed in O- 
type stars of all luminosity classes JChlebowski & Garm anv 
Il99lh . To reproduce the O vi and the N v lines, we have ac- 
counted for shock-generated X-ray emission in the final mod- 
eling stage. Because the two selected stars have not been ob- 
served with X-ray satellites, we have assumed the luminosity 
ratios, logLx/^bob measured for stars of same spectral sub- 
type and luminosity class. For HP 190429 A, we hav e adopted 
logLx/^boi = -7.1, quoted by Pauld rach et alJ d200 lh for the 



04 I(n)f star £ Puppis, whose spectral and physical properties 
are expected to be quite similar to those of HD 190429 A; for 
HD 96715, we have adopted log Ly /Lhni = -6.6, measured f or 
the 04V((f)) star HD 46223 lchlebowsk f& Garmanvll99ll) . 

4. Spectral analysis 

Following the methodology outlined in iBouret et alJ J2003), 
we determined initial estimates of the stellar parameter s 
with Tlusty model atmospheres jLanz & Hubenvl 120031) . 
Subsequent analysis is performed with CMFGEN, using 
Tlusty parameters as initial inputs. At this stage, we still 
allowed for changes in the photospheric parameters, taking 
advantage of the realistic description of photospheric layers 
by CMFGEN. We always checked for consistency between 
Tlusty and CMFGEN determinations. 

Hydrogen and helium lines in the optical spectrum are the 
classical diagnoses used in spectral analyses of O stars. In this 
paper, we presen t an analysis based on the far-ultraviolet spec- 
trum exclusivelv. lHeap e t al. (2005) discussed in detail the UV 
lines and ionization equilibria that are good indicators of stel- 
lar parameters (effective temperature, surface gravity). Table [3] 
lists the most important lines used in our analysis and the quan- 
tities they are most sensitive to. 

As a reference , we a dopted the solar abundances from 
iGrevesse & Sauval Jl998l) . However, we note that significant 
downward revisions of the abundance of light elements in 
the solar photosphere have been recently proposed on the 
basis o f a 3-D hydrodyna mical model of the solar atmo- 
sphere ( Asplu nd et alJl2004T) . bringing the latter values in bet- 
ter agreement with surface abundances of B-type stars in the 
solar neighborhood. In this paper, all chemical abundances are 
quoted by number dens ity relative to hydrogen, or relative to 
IGrevesse & Sauvall Jl998l) solar values. 

4.1. HD 1 90429 A- 04 If + 

The IUE and FUSE spectra offer a broad range of ioniza- 
tion stages and diversity of species that may be used to con- 
strain the effective temperature: He n, C m, C iv, N m, N ry , 

5 iv, S v, S vi, Fe iv, Fe v, Fe vi. Following Heap etaD 12005), 
we used in priority ratios between successive ions of C, N, 
O and Fe. Both the C mAl 176/C wAl 169 and the Fe iv/Fe v 1 
line ratios indicate T e g£ 40, 000 K. The Fe vi lines in the 1250- 
1350 A range are quite weak at these temperatures. Moreover, 
the N 111/1/198 1-991 lines are found to be very sensitive to 
J e ff, as noticed by ICrowfher et alJ ( l2002bl) . and they set a 
lower limit T e g ^ 37, 500 K. The ionization ratio, O iv/l/11338- 
1343/0 v/11371, further supports this lower limit on r e ff. On 
the other hand, Si rv/LU 393- 1402, which are also very sensitive 
to temperature, suggest T e n as 35, 000 K. It is unlikely, how- 
ever, that a contamination of the IUE spectrum by the com- 
panion is responsible for this lower value. Indeed, the large 
IUE aperture contains both the primary and the companion 
HD 190429B (at 2"), but the 09.5 II spectral type of the com- 

1 Fe iv lines between 1500 and 1700 A, and Fe v lines between 1300 
and 1500 A. 
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Table 3. Main spectral lines used in the analysis of the photo- 
spheric and wind properties of HD 190429A and HD 96715, 
and their sensitivity to stellar parameters. 



Dependencies 
Line A (A) HD 190429 A HD 96715 



Siv 


1073-99 


r eff , M, S/H 




P V 


1118-28 


M, /, P/H 


ft. P/H 


m 


1150-54 


ft, O/H 


ft,0/H 


C iv 


1169 


r cff , M, C/H 


r eff , ft, c/h 


C m 


1176 


r cff , M, C/H 


r eff , ft, c/h 


N m 


1182-84 


ft, N/H 


ft, N/H 


Fe vi 


1250-1350 


r cff , ft, Fe/H 


r cff , ft, Fe/H 


O iv 


1338-43 


r cff , M, O/H 


r cff , M, O/H 


O v 


1371 


r eff , M, /, O/H 


r eff , M, /, O/H 


Fe v 


1350-1500 


T eff , ft, Fe/H 


T eff , ft, Fe/H 


Si iv 


1393,1402 


r eff , m, / 




Cm 


1426-28 




C/H, Jeff 


S v 


1502 


M, f , S/H 


ft. 7W 


Civ 


1548-50 


v«,,/5,ft 


v„,Aft 


Fe iv 


1500-1750 


T eff , ft, Fe/H 


T eff , ft, Fe/H 


Hen 


1640 


M, /, He/H 


M, He/H 


Niv 


1718 


M, /, N/H 


M, /, N/H 


Nm 


1748-52 


ft, N/H 


ft, N/H 


Civ 


1860 


C/H, Teff 




Cm 


1875-78 


c/h, r eff 




Nm 


1885 


ft, N/H 





pamon ( Walb orn & Howarfhl2000l) implies that its contribution 
remains small (the luminosity ratio should be about 5). Overall, 
we found that the best fit is achieved for T e ff — 39,000 K, with 
an uncertainty better than 1000 K. This value is in exc ellent 
agreement with the results from iMarkova et alJ d2004l) . who 
used a calibration of T e ff with the spectral type based on re- 
sults obtained from the a nalysis of hydrogen a nd helium lines 
in th e optical spectrum dRepolust et alJl2004l) . On the other 
hand. lGarcia & Bianchil d2004l) derived a slightly lower value, 
T e ff= 37,500 K, based on UV lines formed in the wind. These 
lines are sensitive to the detail of the wind ionization, hence 
to wind clumping that Garcia & Bianchi neglected. Moreover, 
they used carbon and oxygen lines assuming solar abundances 
while admitting that these two species might have highly non- 
solar abundances (see below). 

Our initial estimate of the surface gravity comes from the 
relation between th e spectral type and log g for luminosity 
class I stars (Fig. 2. IMarkova et alJl2004l) . This relation yields 
logg = 3.65 for an 04 supergiant. We then explored models 
with logg ranging from 3.5 to 3.75 (by steps of 0.05 dex, for 
r eff fixed at 40,000 K). We found that the shape of the ob- 
served spectral energy distribution is indeed best reproduced 
for log g * 3.6 - 3.65 (see lLanz & Hubenvll2003l for the de- 
pendence of the SED as a function of log g, at a given r e ff). We 
have adopted, \ogg - 3.6, for the rest of the analysis. 

The visual photometry listed in Table ^ pertains to both 
components A+B. Based on the Hip parcos magni t ude di ffer- 
ence of the two components, Walborn & Howarth (2000) ob- 
tained V = 7.12 for component A alone. The absolute visual 
magnitude is then derived, ass uming a distance d = 2.29 kpc 
for the CygOB3 association (Humphreys 1978). The stellar 



luminosity is calculated by applying a bolometric correction 
from Lanz & Hubenv (2003). The derivation of the stellar ra- 
dius then follows straightforwardly. 

The microturbulent velocity, ft, is determined from the 
iron line strengths (Fe iv and Fe v). The iron abundance is 
kept to the s olar v alue through out this a nalysi s. As noted in 
iHillier et all d2003h and in bouret et alJ fe003l) . O rv/U1338- 
1343 is also sensitive to ft. However, the oxygen abundance 
may have a non-solar value at the surface of an 04 If+ star 
(see below). A consistent fit to the aforementioned lines is ob- 
tained for ft = 15kms _1 . We note that these line s have been 
simila rly used in other studies dBouret et alJl2003l iHeap et all 
2005), yielding a range of microturbulence values for O stars 
(ft = 2-25 km s~'). It is therefore unlikely that the high micro- 
turbulence derived here results from inadequate model atoms. 

The helium abundance is poorly constrained by the FUV 
spectrum alone. We first used y - He/H = 0.1 (by number), 
but subsequently we increased the abundance to y = 0.2. The 
higher helium abundance provides a better fit to He n/11640 
(once the mass loss rate is fixed at the value derived from 
other lines). The photospheric He H/U085 line, which is not 
affected by the wind in our models, could not be used be- 
cause of strong blends from IS N n lines. The enhanced helium 
abundance is consistent with the evolutionary status of an 04 
supergiant dMevnet & Maede3l2000l) . but we cannot rule out 
that it may be further revised when optical data is analyzed in 
addition to the FUV spectrum. A simultaneous fit of all He i 
and He n lines througho ut the spectrum is, however, rarely (if 
ever) achieved (see, e.g.. lHillier et alJl2003l) . Additionally, we 
note that changes in the helium abundance impact the predicted 
P v/Ul 1 18-1 128 lines, because the ionization threshold of P v 
is at 191 A, close to the He n Lyman limit. Models with y = 0.2 
yield a phosphorus abundance that is closer to the solar value 
compared to values derived with y — 0.1 models. On the basis 
of these various arguments, we finally adopted y = 0.2. 

The carbon abundance is constrained by the Ciii/ill76, 
C ni/11923, and Civ/11169 photospheric lines. Only models 
strongly depleted in carbon match these lines. A low car- 
bon abundance is also consistent with the C inA911 and the 
C rv/l/U548-1550 resonance lines, although these lines are not 
sensitive indicators because they show saturated P Cygni pro- 
files that are only weakly sensitive to the wind velocity law. 
We also checked that the models do not predict a strong 
C m/1 4647-5 1 emission that is no t observed in the optical spec- 
trum (Walborn & Howarth 2000). The best fit to the lines dis- 
cussed above is obtained for C/C Q = 0.05. This value is consis- 
tent with the evolved nature of HD 190429 A. 

All nitrogen FUV lines, but N iv/11718, show that the sur- 
face of HD 190429 A has been enriched in nitrogen. This con- 
clusion is additionally su pported by the observed N mAA4634- 
4640 emission (see, e. g.. lWaiborn & HowarthfcO OO). We have 
derived a nitrogen surface abundance, N/N Q = 4.0. The ni- 
trogen enrichment of HD 190429 A is consistent with its be- 
ing highly evolved, as significant nitrogen enrichment together 
with carbon depletion is predicted by stellar evolution mod- 
els for stars in the hel ium burning ph ase. Note in this respect 
fhat lConfietaHdl995l) and lWalborn & Howarfhld2000l) already 
suggested an advanced evolutionary status for HD 190429 A 
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Fig. 1. Best fit to the photospheric lines used to derive the stellar parameters of HD 190429A (T e ff- 39,000 K, logg = 3.6, 
15 km s _1 ). The adopted abundances are listed in Tabled 



based on K-band and optical data (respectively), which they 
interpreted as indications that the star is well advanced toward 
the WN stage. 

The lack of wind-insensitive oxygen lines hampers an ac- 
curate determination of the oxygen abundance. The strong 
O vi/l/11032-1038 are fully controlled by the X-ray flux and 
other wind parameters. The O V/11371 line is rather weak and 
is most sensitive to the mass loss rate and clumping parame- 
ters. Finally, although mostly formed in the photosphere, the 
O ivAA 1338- 1343 lines show a blue asymmetry that indicates 
a contribution from the wind. Assuming g t - 15kms~' (see 
above), we had to adopt a low oxygen abundance, O/O = 
0.1, in order to match the O IV/LU338-1343 lines. However, 
we always failed to reproduce the O v/11371 line with homoge- 
neous wind models. This low abundance is also supported by 
the weak O m triplet (/HI 150-1 154 A) in the FUSE spectrum. 

Phosphorus and sulfur abundances are not expected to be 
affected by nucleosynthetic processes. We thus adopted ini- 
tially (and kept for most of the analysis) solar abundances. 
There are, however, indications that both elements might 
be slightly depleted. For homogeneous wind models, a low 
phosphorus abundance, P/P = 0.1, is required to match the 
P vAAl 1 18-1 128 resonance doublet. On the other hand, we de- 
rive P/P = 0.5 from a clumped model with a filling factor 



fa, = 0.04 (see below). All sulfur ions indicate that S/S Q = 0.5 
for homogeneous mode ls and S/S n = 0.9 f or clumped models 
(/ M = 0.04). Similarly. Ikillier et all ^2003) derived lower than 
one-fifth solar-scaled P and S abundances for the SMC 07 Iaf+ 
supergiant AV 83. However, given the r ecent downward revi - 
sion of the solar P and S abundances ( Asplu nd et al.lE004l) . 
we argue that the abundances derived from the clumped wind 
model roughly remain in agreement with solar values within 
the uncertainties (0.1 to 0.2 dex). 

A detailed examination of the complete FUV spectrum re- 
veals that only a few lines are influenced by the stellar wind, 
even for a supergiant with such an early spectral type and 
strong wind (M as few 10~ 6 M o yr -1 ). In particular, these lines 
include a few strong Fe lines whose cores are filled a lit- 
tle by wind emission. For most other lines, we have a very 
good consistency between Tlusty and CMFGEN model spec- 
tra. This agreement demonstrates that we may use NLTE pho- 
tospheric models to derive reliable stellar parameters and abun- 
dances, even for Galactic O supergiants with relatively dense 
winds. Furthermore, we stress that we had to use large (£ f = 
15kms _I ) microturbulent velocities in CMFGEN wind mod- 
els too. Despite an early comment by iKudrit/kij |l992|, this 
demonstrates that the need for microturbulence is not restricted 
to hydrostatic models that neglect the atmospheric velocity 
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Fig. 2. Best fit to HD 190429 A wind-sensitive lines, obtained with clumped (full grey line) and homogeneous (dashed grey line) 
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. The short wavelength side of N V/U240 is affected by interstellar Lya absorption (see Table 1 for n(H i)). Notice in particular the excellent fit 

to both lines of the P v doublet achieved with the clumped wind model. 



field. We display in Fig^the best CMFGEN model fit to some 
of the lines discussed above, which are only weakly sensitive 
to the wind properties. 

We now turn to the analysis of the wind properties of 
HD 190429 A. The parameters of the /6-type velocity law are 
derived from the saturated C iv/l/11548-1550 doublet and from 
the N iv/11718 P Cygni profile. A good match is achieved for 
j6 = 0.8 and v ro = 2300 km s . We have used an outward in- 
creasing microturbulence velocity from the photospheric value 
(15kms _1 ) to a maximum of 200 kms -1 . 

To derive the mass loss rate, we started w ith an estimate 
based on Her modeling, M = 1.4 10~ 5 M Q yr~' (Markova et al. 
120041) . Such a high value can be excluded from the FUV 
spectrum, because a wind model assuming this value pre- 
dicts much stronger P Cygni profiles than those observed. 
Furthermore, this model also produces asymmetric Fe iv and 
Fe v lines with a b l ue ex tended absorption that is not observed. 
iGarcia & Bianchil J2004I) derived such a high mass loss rate 
value from the UV wind lines. However, from their Fig. 7, we 
see that their model predicts C mvU176 and O iv/11340 emis- 
sions that are not observed. Generally, all other wind lines are 



also predicted too strong by their model. We then built a series 
of wind models, decreasing M in several steps. Using homoge- 
neous wind models, we have not been able to achieve a good 
match to most observed P Cygni profiles and to other wind sen- 
sitive lines like O IV/L11338-1343 and He ilH640. 

Therefore, we switched to clumped wind models with 
the p resumptio n tha t the high mass loss rate derived by 
Mar kova et alJ J2004) was the result of the ef fect of clumps 
on Ha line formation (as already suggested by Rep olust et alJ 
12004) . We searched thus for a combination of M and /«,, keep- 
ing M/yffZ ~ 1.4 lO^Moyr 1 . The best match to the FUV 
lines was obtained for M = 1.8 10 _6 M o yr _1 and f m = 0.04 
(see Fig. The fit to P v/Ul 1 18-1 128 and to O v^1371 line 
profiles are the most striking improvements achieved with the 
clumped wind model. The agreement with the observed pro- 
file gets also better for N iv/11718, while there are little pro- 
file changes for the C iv resonance doublet or He H/11640 com- 
pared to profiles predicted by a homogeneous wind model. The 
sensitivity of P v lines to clum ping was first e stablished by 
ICrowther et ail J2002bl) and iHillier etall J2003I) . We defer to 
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Fig. 3. Best fit to photospheric lines used to derive HD 96715 stellar parameters (r e fj = 43,500 K, logg = 4.0, = 15 kms '). The 
adopted abundances are listed in Tabled 



13 a discussion of the physical implications of clumping on 
line formation. 

Other wind line profiles, such as O vi/11036 and N V/11240 
are strongly influenced by X-ray wind emission (see 35.21 and 
Fig. Q, and provide little additional constraints on the mass loss 
rate or on the clumping properties of the wind of HD 190429 A. 
P Cygni profiles of lower ions, like C ni/1977 or N ni/1991, are 
severely blended with H2 lines that hamper an accurate deter- 
mination of M from these lines. 



4.2. HD 96715 - 04 V((f)) 

The IUE spectrum of HD 96715 exhibits conspicuous 
lines of Fevi (1250-1350 A), Fe v (1300-1550 A) and Ferv 
(1550-1700 A), as well as O IV/U1338-1343, C mXU 426- 1428 
and S V/11502. The FUSE spectrum shows Civ/11169 and 
Cm/11176 lines, as well as O hlUI 150-1 154, S IV/U1073- 
1099. The lines are relatively narrow, implying an apparent ro- 
tational velocity, v sin i = 80 km s _1 , and resulting in less severe 
blending problems from iron lines. 

Using the ionization balance between these ions, we de- 
rived a best overall match with a Tlusty model having T e s 
= 43,500 K and logg = 4.0. A realistic estimate of the uncer- 
tainty on the r e ff determination is about ± 1500 K. On the other 



hand, the sensitivity of the FUV spectrum to logg is small. We 
can, however, firmly exclude models with gravities differing 
by ± 0.25 dex in log g, based on the comparison of theoretical 
spectral energy distributions with the IUE spectrum corrected 
for reddening (assuming E(B - V) = 0.42, Table [Q. The best 
overall agreement is obtained for models with log g = 3.9 + 0. 1 . 
The relation between s pectral type and sur face gravity for lu- 
minosity class V stars (Mar kova et aljE004l) yields logg= 3.9 
in agreement with our estimate within the error bars. We have 
finally adopted logg= 4.0, see Table@] 



iGarcia & Bianchil ll2004) obtained a markedly lower value, 
T e ff= 39,000+2,000 K. They established upper and lower tem- 
perature limits from the absence of Pv/11118, 1128 and 
Ov/11371. In particular, the weakness of the Ov line con- 
strains r e ff<40,000K. Homogeneous wind models predict a 
very strong O v feature at temperatures higher than 40 kK, thus 
forcing an artificially low T e s to be adopted. Garcia & Bianchi 
methodology yields temperatures that are systematically lower 
than all other recent analys es o f optical and UV spe ctra of O 
stars (see lHeap et alJ J2005h and lMartins et alJ (|?005) who pro- 
vide a detailed comparison of these studies and who are led 
to exclude Garcia & Bianchi results). We believe that their 
methodology underestimates temperatures because of the re- 
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liance on wind lines that are sensitive to the wind structure and 
abundances as well as to T e g. 

We adopted the abso lute magnitude, M v = -5.5, from 
iHowarth & Prinial (1989), and derived the stellar luminosity 
and ra dius using a bolometric correction from lLanz & Hu benv 
The stellar mass then follows from logg and R*. 

The microturbulent velocity has been derived from O iv, 
Fe iv, Fe v, Fe vi, and S v lines. The best match is ob- 
tained for = 15 km s" 1 , assuming solar abundances for these 
species, consistently with the expected evolutionary status of 
HD 96715. 

We started the analysis using a normal helium abundance, 
y = 0.1. Subsequently, we slightly decreased it to y = 0.09 to 
improve the match to He II/U085 and He ii/il640. Admittedly, 
this value should be considered as a provisional estimate un- 
til an analysis of the optical spectrum becomes available. Yet, 
there is no evidence of helium enrichment at the surface of 
HD 96715. 

The carbon abundance was derived from the photospheric 
lines, Civ^1169, C mil 176, and C IIIA1426-1428. Using a 
Tlusty model with the stellar parameters previously derived, 
we found that C/C Q = 0.5 matches well these carbon lines. 

Like HD 190429 A we constrained the nitrogen abundance 
from N m/t/tl 182-1 184 and N mXil748-1752. We found that 
an enhancement of a factor of 4 with respect to the solar value is 
required to match the observed lines. Note that wind lines such 
as N iV/11718 further support this determination since no match 
can be achieved without a significant nitrogen enrichment. Our 
result is also consistent with a line blend, N nAA 1084- 1086, al- 
though the presence of superimposed IS lines makes the deter- 
mination of the photospheric contribution difficult. 



To measure the oxygen abundance, we relied on 
O mAAl 150-1 154 photospheric lines, which are clearly seen 
in the FUSE spectrum. A very good fit to this triplet is ob- 
tained for O/O = 0.9, which also results in a good fit of 
O iv/l/11338-1343. Fig. [5] displays the best model fit achieved 
with CMFGEN to the photospheric lines discussed in this sec- 
tion. 

The wind velocity law has been determined on the satu- 
rated P Cygni profile of C iv/l/11548-1550. The terminal veloc- 
ity and maximum turbulent velocity have been measured on 
the blue-side of the absorption component; we have adopted 
a terminal wind velocity, v m = 3000 km s -1 , and a maximum 
turbulent velocity in the wind of 250 km s . The P Cygni pro- 
file of C IV/U1548-1550 was preferred to N V/11240 because 
the latter is affected by a significant contamination due to IS 
Lya and because of its well-known sensitivity to X-rays. The 
wind acceleration parameter, /?, and the mass loss rate have 
been inferred from the P Cygni profiles of C iv/l/11548-1550 
and N iv A 17 18 together with other wind sensitive lines like 
O ivi/11338-1343 and He nA1640. For the adopted abundances, 
we found that a homogeneous wind model with f3 = 0.8 and M 
= 1.8 10~ 6 M o yr~' reproduces well these lines, although the 
match to N IV/U718 remains poor. 

On the other hand, the O v/11371 line profile predicted by 
homogeneous wind models is stronger than observed. This be- 
havior has been found in all previous studies of O-type stars 
that show O V/H371. More specifically, the wind models pre- 
dict too strong absorption at moderate velocities, too little ab- 
sorption at low velocities, and a too strong redward emis- 
sion. Homogeneous wind models thus fail to reproduce this 
line a s in our earlier st udy of main-sequence O stars in the 
SMC JBouret e t al. 2003). Therefore, we have calculated sev- 
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eral models with different clump volume filling factors /«,, ad- 
justing the mass loss rate so as to maintain a good fit to C rv 
and N iv P Cygni profiles. Fig.0]shows that we can achieve a 
good match of O V/U371 when assuming a very small filling 
factor, f x = 0.02, with a mass loss rate, M- 2.5 10~ 7 M yr -1 . 
A larger filling factor, /«, = 0.1, and a lower oxygen abun- 
dance by a factor of 4 also improves over the case of a homo- 
geneous wind model, but does not fit the observed profile as 
well as a model with the smaller filling factor. Larger volume 
filling factors may seem more reasonable or, at least, they are 
more in line wi th earlier studies of Wolf-Rayet clumped winds 
faamann & Koesterkel 19981) and of SMC O stars feouret et alJ 
2003|. However, other oxygen lines, like O m/Lll 150-1 154 and 
O iv/l 1338- 1343, are no longer fitted with the lower oxygen 
abundance. These lines are indeed sensitive to changes in the 
adopted oxygen abundance since they are unsaturated. There is 
no compelling reason for adopting an oxygen abundance lower 
by a factor of 4 relative to the solar value, especially since car- 
bon is not depleted as well. Carbon would indeed be expected 
to show first a depletion typical of CNO-cycle processed ma- 
terial. We have therefore adopted the clumped model with the 
very small volume filling factor. 

The steep transition between the absorption and emission 
components indicates that clumps start forming at low veloci- 
ties, v c i » 30kms~', just above the sonic point. At this point, 
we discovered that N iv/11718 behaves similarly to O V/11371, 
and is thus sensitive to wind clumping too. The fit to N IV/11718 
is improved by using clumped wind models, especially at the 
transition between the absorption and emission components. 
The extended blue absorption is affected similarly to O vA 1 37 1 . 
The nitrogen abundance has been kept to the value derived from 
photospheric lines. Finally, we stress that we have achieved a 
good match of both the N iv and the O v lines using the same 
low filling factor for the clumps, which thus provides convinc- 
ing support for the highly clumped model. The profile of other 
lines (C iv and N v) is not directly affected by clumping. 

5. Discussion 

Table |4] summarizes the results of our analysis. Conservative 
estimates of the uncertainties are +5% on T e ff, +0.1 dex on 
logg, +2 to iSkms 1 on£ (for HD 190429AandHD 96715, 
respectively), and +0.1 to ±0.2 dex on chemical abundances. 
Abundances are listed as number densities and refer to values 
adopted in clumped wind models (lower abundances might be 
preferred with homogeneous wind models in a few instances, 
see §4). The mass loss rates are derived with an accuracy of 
10-20 % from the fitting process, though this does not include 
systematic errors. We discuss now the implications of our re- 
sults, and compare them to other observational and theoretical 
studies. 

5.1. Evidence of clumping from the FUV spectrum 

^From the foregoing analysis, we conclude that Pv/1/11118- 
1128 and Ov/11371 provide the best diagnostics of clump- 
ing in the wind of O s t ars, supporting the s emin al studies of 
ICrowther etail J2002bl) . iHillier et all J2003h and bouret et ail 



Table 4. Stellar parameters, wind properties, and chemical 
abundances for HD 190429 A and HD 96715. The listed abun- 
dances correspond to values derived from clumped models (see 
§4.1). 



Star 


HD 190429 A 


HD 96715 


Spectral Type 


04If+ 


04 V((f)) 


7W[K] 


39000 


43500 


logg (cgs) 


3.6 


4.0 


L[L G ] 


7.9 x 10 5 


4.6 x 10 5 


ft [kms- 1 ] 


15 


15 


v sin i [km s _1 ] 


160 


80 


Homogeneous winds 


M [M Q yr'] 


6.0 x 10" 6 


1.8 x 10 -6 


Voo [kms 1 ] 


2300 


3000 


P 


0.8 


1.0 


Clumped winds 


M [M yr>] 


1.8 x 10- 6 


2.5 x 10- ' 




0.04 


0.02 


Surface abundances 3 


.V (He/H) 


0.2 


0.09 


c/c Q 


0.05 


0.5 


N/N Q 


4.0 


4.0 


o/o Q 


0.1 


0.9 


Si/Si Q 


1.0 


1.0 


P/Pe 


0.5 


1.0 


s/s Q 


0.9 


1.0 


Fe/Fe G 


1.0 


1.0 


N/C b 


80 


8 



a Solar abundances from lGrevesse & Sauva 

b Abundance ratio relative to the solar ratio, N/C = 0.25. 

ll2003t) . and extending them to the case of Galactic stars. In at- 
tempting to reproduce the O v/11371 line profile, Bouret et al. 
explored a number of alternative explanations besides clump- 
ing. They emphasized that the key point consists in predicting 
the correct O v/O iv ionization structure in the vicinity of the 
sonic point. To alter ionization, they considered X-ray emission 
from shocks in the wind (Lx/L\, \ = 2 10~ 7 ), adiabatic c ooling, 
larger model atoms, and alternate recombination rates ( Nahar 
1999). All these changes resulted in essentially no difference in 
the model spectrum, and wind clumping remained as the only 
viable explanation. 

The case of P v is different. ICrowther et alJ d2002bl) and 
HillieretalJ ll2003l) showed that an effect on the line profile 
similar to clumping could be obtained by drastically reducing 
the phosphorus abundance because P v reaches its maximum 
and is the dominant ion th roughout the wind of late O super- 
giants. Hillieret al. (2003) achieved a good match to the ob- 
served P v doublet in the SMC 07Iaf+ supergiant AV 83 with 
homogeneous wind models using a low abundance P/P = 0.04 
(compared to the adopted overall SMC metallicity, Z/Z = 
0.2). While there are remaining uncertainties relative to the 
phosphorus abundance baseline and the chemical evolution 
of the Magellanic Clouds leaving th us some leeway toward 
a definitive conclusion, we note that Pauldrach et a U (11994 
2001) similarly claimed that two Galactic supergiants, namely 
( Puppis (04 I(n)f) and a Cam (09.5 la) must have a very low 
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Fig. 5. Phosphorus ionization fractions in clumped (black lines) 
and homogeneous (grey lines) wind models of HD 190429 A. 
The P iv fraction remains smaller than 1 % throughout the wind 
(max. at v(r) » lkms~'). Notice that the velocity scale is 
shifted toward higher velocities for the homogeneous wind 
model because of a higher mass loss rate (the maximum of 
P v occurs at similar densities in the two models), and that the 
adopted phosphorus abundances are different in the two models 
(P/P = 0.1 and 0.5, for the homogeneous and clumped models, 
see §4.1). 

phosphorus abundance based on their analysis of Copernicus 
data. In the case of Galactic stars, a low phosphorus abundance 
seems a much less likely explanation of the weak P v resonance 
lines, and we thus speculate that Pauldrach et al.'s results actu- 
ally indicate that ( Puppis and a Cam also possess a highly- 
clumped wind. Our study shows that clumping combined with 
a slightly reduced phosphorus abundance is required to match 
the P v resonance doublet in HD 190429 A, wh ile still fitting 
well the rest of the spectrum. Hill ier et al J (12003 ) also achieved 
an excellent match of the full spectrum (from FUV to the opti- 
cal) of AV 83 using / M = 0.1 and P/P = 0.08, but the evidence 
for a lower phosphorus abundance remains relatively weak be- 
cause of the possibility of adopting a lower clump filling factor 
coupled to a higher abundance. These results, therefore, only 
hint at a depletion of P in O supergiants with respect to the 
other metal abundances. 

In addition to O v and P v, we finally found that N iV/il718 
is also an indicator of wind clumping. The N iv line profile 
exhibits the same behavior as Ov/11371. The two lines cor- 
respond to the same transition, 2s( 2 S)2p 'P - 2p 2 1 D, in an 
isoelectronic sequence. The O v and the N iv lines can be fitted 
with a clumped wind model having the same clumping param- 
eters. 

For both stars, we found small volume filling factors, = 
0.04 and 0.02, for HD 190429A and HD 96715, respectively. 
Although these small values may seem problematical, they 
are supported by other recent results ab out wind clumping . 
In the analysis of 25 O stars in the LMC, [Massa et al. (2003) 
found that the P v ionization fraction never exceeds 0.20, al- 
though it is expected to be the dominant stage in some of 



Fig. 6. Oxygen ionization fractions in clumped (black lines) 
and homogeneous (grey lines) wind models of HD 96715. 
Notice that the velocity scale is shifted toward higher veloci- 
ties for the homogeneous wind model because of a higher mass 
loss rate. 

the stars. They argued that this implies either that the cal- 
culated mass loss rates or the phosphorus abundance are too 
large, or that the winds are strongly clumped (J^ < 0.05). 
The first evidences of wind clumping were found in WR 
stars for which small filling factors (f ro ~ 0.1) were de- 
rived too jHamann & Koesterkdl 19981 .[Hillier & Millerlll999l. 
ICrowther et alJ2002al) . lKurosawa etail (l2002h and lSchild etall 
(2004) derived even lower values, * 0.05 - 0.075, for the 
WR stars V444 Cyg and y 2 Vel. Hydrodynamical 1-D calcu- 
lations of line-driven wind instabilities support low volume 
filling factors (/„ < 0. 1 : iRunacres & Owockil l2002). but ini- 
tial results of a restricted 2-D app roach indicate larger values 
(/«, * 0.2:lDessart & Owockil2003l). 

5.2. Effect of clumping on wind ionization 

We now examine the reasons why these line profiles are mod- 
ified by clumping. Essentially, the presence of optically thick 
clumps separated by transparent voids has two major conse- 
quences: (i) the wind over-densities strengthen the emission of 
density-sensitive lines; (ii) the number of recombinations in- 
creases because of the higher density in the clumps, thus re- 
ducing ionization. The competition between these two effects 
may either increase or decrease emission. Because of the ra- 
dial dependence of the volume filling factor, lines of different 
strength behave in different ways. This is shown in Fig.|2]and 
Fig.|4]where profiles for clumped and homogeneous models are 
compared for HD 190429 A and HD 96715. 

In 07 supergiants, P v is the dominant ionizati on stage of 
phosphorus throughout the wind jHillier et al.l2 003). However, 
HD 190429 A being a hotter supergiant, P v (in the photo- 
sphere) and P vi (in the wind) are now the main ionization 
stages. The P iv fraction always remains small (less than 1%). 
Figure |5]illustrates the effect of clumps in the stellar wind: we 
find a lower ionization, P vi/P v, compared to the homogeneous 
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wind model, demonstrating the increased recombination from 
P vi to P v in the clumps. In HD 190429A, the P vAAl 1 18-1 128 
resonance lines reflect the factor of 2 between the /-values of 
the two lines, thus indicating that these lines are unsaturated. 
The interpretation of the strength of these lines is, however, not 
straightforward, that is, it depends not only on the phospho- 
rus density (a product of the mass loss rate and the phosphorus 
abundance), but it also depends on the detail of the wind ioniza- 
tion. The dependence with the clumping parameter, /«,, is in- 
direct through the requirement that the ratio Mj JfZ has to re- 
main constant for keeping a good match to other lines, such as 
Ha and C iv/il550. A low phosphorus abundance (P/P Q = 0.1) 
must be adopted with a homogeneous wind model, thus imply- 
ing significant clumping (low /«,, hence low M) is necessary 
to increase t he ph osphorus abundance to a near-solar value. 
iMassa et alJ J2003) followed the same argument to explain the 
systematically low P v ionization fractions (<20%) that they 
derived for 25 O stars in the LMC. We stress that our clumped 
wind model matches the two P v resonance lines very well, pro- 
viding a strong evidence that the wind ionization is correct and 
that there is no particular issue such as covering factors with 
our treatment of clumps. 

The oxygen lines, however, respond in a different way to 
clumping. They are affected by clumping via a lower ionization 
in the wind - see Fig. [6] for HD 96715. Notice that ionization 
in fact is not lower in the photosphere as Fig.[6]would suggest: 
there is a shift in velocity space because of the higher mass loss 
rate of the homogeneous wind model. At similar photospheric 
densities, the ionization is the same. On the other hand, O iv 
is the dominant ion throughout the wind of HD 190429 A. The 
homogeneous model predicts an increase of O v at moderate 
velocities (over 30% atfew 100'skms -1 ), but increased recom- 
bination due to the presence of clumps keeps the O v fraction 
below a few percent in the clumped wind. The lower O v frac- 
tion in the clumped wind models weakens the blue-shifted line 
absorption substantially. The wind ionization is sensitive to the 
clumping parameter, /«,, and O V/11371 is the best indicator of 
clum ping in the wind of HD 96715, similarly to the SMC 04 
stars jBouret et alJ2003t) . 

Besides the O v and Pv lines, we found that N iV/11718 
is also sensitive to clumping. It is a special interest because 
clumping decreases N v/N iv ionization in a way quite similar 
to O v/O iv. N iv is the dominant ion throughout the wind of 
HD 190429 A, while it is the dominant ion only in the photo- 
sphere of HD 96715 (N v is the dominant stage in the wind). 
The N iv/11718 responds to clumping like O V/11371. The ab- 
sorption component of NiV/11718 is globally weaker in the 
clumped model, with the exception of velocities between « 
15kms _1 and 650km s _1 (Fig.|4). The direct mapping of the 
velocity law on the absorption component implies that these 
velocities correspond to the region where N v/N iv ionization 
is lower in the clumped model. In the outer wind, at higher 
velocities, N v recombines into N iv in the homogeneous and 
clumped models of HD 96715 because of lower temperatures. 
The predicted profile at these high velocities (v > 650km s _1 ) 
then reflects the higher mass loss rate adopted for the homoge- 
neous wind model. 




1230 1240 1250 

Wavelength [A] 



Fig. 7. Effect of X-ray wind emission on the predicted 
N V/11240 line profile of clumped wind models of HD 190429 A 
(full line: logLx/^boi = -7.1; dotted line: no X-rays). All the 
other lines sensitive to wind clumping and shown in Fig.|2]are 
not affected by X-ray emission. 



Several other lines are displayed in Fig. 13 showing smaller 
changes because of clumping. For instance, S V/11502 shows a 
slightly stronger blueshifted wing because of increased recom- 
bination from S vi to S v. On the other hand, the C iv resonance 
doublet is only marginally sensitive to clumping because these 
lines are saturated. 

X-ray emission in the wind also changes the wind ioniza- 
tion structure, increasing the populations of "superions" like 
N v, O vi, and S vi. We found that the profiles of N V/11240 and 

vi/11036 are substantially changed by X-ray emission (see 
Fig- EI m e O vi lines are not shown because of strong interstel- 
lar H2 absorption). All other FUV lines shown in Fig. 13 espe- 
cially those sensitive to wind clumping, remain unaffected. In 
the outer wind of HD 190429A (v(r) > 500 km s 1 ), the N v 
is increased by a factor of several up to several orders of mag- 
nitude at terminal velocity by the X-ray wind-shock emission 
compared to predictions neglecting X-rays. In the inner wind 
(v(r) < 500kms _1 ), N v increase is more modest (+20% com- 
pared to predictions without X-rays). The N v fraction reaches 
about 10% in the inner wind and always remains larger than 

1 % in the outer wind while it precipitously drops in the model 
without X-rays. This change in the wind ionization structure 
results in the marked line profile change displayed in Fig.0 

Finally, we need to stress here that clumping is not a sim- 
ple numerical trick introduced to correct the ionization bal- 
ance of a single species (e.g., O v/O iv). We showed here that 
our clumped models change in a consistent way the ioniza- 
tion balance of different species, with different ionization po- 
tentials, and affect the corresponding lines differently as they 
form in different parts of the wind and respond differently to 
the nature of clumping in the wind. We therefore conclude that 
our description of wind clumping, albeit simple, must be basi- 
cally correct because of our ability to fit lines of different ions 
and species using the same clumping parameters. The clumped 
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Fig. 8. Predicted Ha profiles for HD 190429 A, calculated with 
clumped (grey line) and homogeneous (black line) wind mod- 
els. The corresponding parameters are listed in Table 4. 



Fig. 9. Predicted Ho- profiles for HD 96715, calculated with 
clumped (grey line) and homogeneous (black line) wind mod- 
els. The corresponding parameters are listed in Table 4. 



wind models therefore provide a sound basis for determining 
improved mass loss rates. 

5.3. Clumping, Ha, and the mass loss rate 

We discuss now the implication of wind clumping on the de- 
rived mass loss rates, comparing our results based on the anal- 
ysis of the FUV spectrum with determinations based on Ha 
and thermal radio emission. Such a comparison is essential be- 
cause the quantity that is actually measured from the FUV lines 
is the product of the mass loss rate by the appropriate ionization 
fraction, and we just argued that clumping results in changing 
the mean ionization in the wind. We believe that CMFGEN 
models provide a reliable description of the ionization struc- 
ture of O star winds thanks to the detailed, unified treatment of 
NLTE metal line blanketing in the photosphere and in the wind. 
Yet, a comparison with results derived from Ha is crucial since 
the latter are much less sensitive to many factors influencing 
the wind ionization. Moreover, the FUV lines behave linearly 
with the wind density while Ha and the radio emission have a 
quadratic dependence. 

Ha- is a recombination line, and is sensitive to the densest 
regions of the wind. Therefore, this line is expected to con- 
strain clumping deep in the wind. On the other hand, thermal 
emission in the radio p robes the outer regions of the wind. 
iKudritzki & Pulsl J2000) argued that Ha should not be sensitive 
to wind clumping contrary to the submillimeter/radio emission 
because they expected clumps to form far out in the wind. This 
is contrary to our findings that clumps already form close to the 
sonic point. 

Although we have no Ha observations in hand, we can 
compare our predictions for homogeneous and clumped wind 
models with published results for HD 190429 A. We display 
in Fig.|S]the predicted Ha profiles calculated with the param- 
eters derived from our analysis (see Table |4). We found only 
little changes between the predicted profile for the homoge- 
neous and clumped wind models, thus indicating that Ha ob- 
servations alone would not exclude one or the other model (but 



different mass loss rates would be inferred). Our predicted line 
profiles are consistent with Walborn & Howarth's (2000; their 
Fig. 2) published profile. 

Table |5] lists the mass loss rates derived from Ha and ra- 
dio observations and compares them to our results. The mass 
loss rates (Ha and radio) depend on the adopted distance, 
proportionally to d 3 ^ 2 (e.g., Scuderi et al. 1998). The val- 
ues listed in Table [5] have been corrected of this dependence. 
The mass loss rate for a homogeneous wind derive d from the 
FUV lines is consistent with ear lier analyses of Ha jLeithererl 
198 8:1: lLamers & Leifhererl 1993b . but lScuderi et alJJl998h and 
Markova et alJ J2004I) derived a value twice as large. Unified 
comoving-frame models were used by Markova et al. and in 
our stud y, while the other an alyses used the Sobolev approxi- 
mation. Scud eri et alJ dl998h also report an increase by a fac- 
tor of 2 of the Ha equivalent width between 1988 and 1991. 
They obtained also a single observation of the free-free emis- 
sion at 8.45 GHz. It is therefore difficult to reach a firm con- 
clusion from this comparison since intrinsic stellar variability 
most likely is at the origin of some of the differences. The main 
conclusions thus remain that (i) the profile from the clumped 
wind model would match the observed Ha profiles as well as 
the predicted profile from the homogeneous wind, and ( ii) the 
derived mass loss rate from the clumped wind model is a factor 
of 3 lower. 

Finally, hydrogen infrared recombination lines such as Bra 
and Bry may be strongly affected by clumping in the wind, be- 
cause of their s ensitivity to the squar e of the local density. As 
pointed out by iLenorzer et alJ (12004), clumping results in in- 
creased recombination rates inside the clumps ((p 2 ) > (p) 2 , as- 
suming an infinite contrast of the clumps with the voids). The 
presence of clumps deep in the wind would thus explain the 
present discrepancies between observed and theoretical profiles 
(and EWs) of lines such as Bry ( ILenorzer et alJl2004 . Indeed, 
clumping close to the photosphere (where Bry forms) would 
increase the collisional excitation rates driving the excitation 
ratios closer to LTE. These recombination lines would thus ex- 
hibit (pure) emission profiles rather than P Cygni profiles (cf. 
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Table 5. Mass loss rates and modified wind momenta measured 
for HD 190429 A and HD 96715. 



Star HD 190429 A HD 96715 



Mass Loss Rates [10 6 M Q yr '] 
This study (far-UV) 



Clumped wind 


1.8+0.3 


0.25±0.05 


Homogeneous wind 


6.0 


1.8 


Ha 






Leitherer (1988) 


9.8 




Lamers & Leitherer (1993) 


5 4 +31 

J - -2.3 




Scuderietal. (1998) 


11.8±1.8 




Markovaetal. (2004) 


14.2 




Radio 




Scuderietal. (1998) 


7.2±0.9 




Modified Wind 


Momentum 




log (M Voo 






Clumped wind 


29.07 


28.21 


Homogeneous wind 


29.66 


29.07 


ReDolust et al. (2004) 


29.54 


29.11 


Vink et al. (2000) 


29.55 


29.13 



IPuls et alii 19961) . An interesting point noted bv lLenorzer et alJ 
is that their predicted Bra agree fairly well with obser- 
vations. This suggests that the degree of clumping might not 
be constant throughout the wind, but rather peaks and then de- 
creases further out in the wind. The location of the strongest 
clumping (e xpected to depend on the m aximum line strength 
according to lRunacres & Owockil J2Q02I) ') might then be in the 
region where Bry forms. Bra which forms over a larger vol- 
ume (because of its larger oscillator strength) would then be 
less affected by clumping. 

Our estimate of M for HD 96715 derived from the homoge- 
neous wind model c ompares well to mass loss rates for stars of 
same spectral type JLamers & Leifhererlll993l: iRepolust et ail 
2004). Ha is always observed in absorption for these stars, 
which is predicted both for the clumped and homogeneous 
wind models. Fig.|9]shows that Ha is somewhat filled by wind 
emission for the homogeneous wind model, but it remains un- 
likely that either model could be excluded based on Ha ob- 
servations. The clumped model yields a mass loss rate that is 7 
times lower than the value derived from the homogeneous wind 
model. 

5.4. The Wind-momentum Luminosity Relation 

The presence of clumps deep in the wind of O stars has been re- 
cently invoked bv lPuls etall J2003I) and Repol ust et al.1 (|2004) 
to explain their findings of two distinct Wind-momentum 
Luminosity Relations (WLR) for stars with Ha in emission 
and in absorption, r espectively . While the WLR for the latter 
objects agrees with IVink et al.l J2000l) theoretical predictions, 
it is clearly separated from the empirical WL R for emission- 
line o bjects by a significan t offset . As shown in Mar kova et alJ 
(2004) and Repolust et al. (2004), this offset practically van- 
ishes when the mass loss rates for the emission-line stars are 
reduced by a factor of 0.44. Because clumped winds mimic 
winds with higher mass loss rates, these authors suggested that 



the higher empirical WLR for emission-line objects reflects the 
influence of clumps in the Ha line formation region. 

From our analysis, we conclude that mass loss rates from 
the homogene ous wind models are very similar to those pre- 
dicted by the IVink et all fcOOOh formula using the parameters 
derived for HD 190429A and HD 96715. Table|5]lists Vink et 
al.'s predicted modified wind momenta together with the val- 
ues derived from our analysis. As might be expected from the 
lower mass loss rates derived from the clumped wind mod- 
els, the wind momenta calculated from these models are sig- 
nificantly low er than the theoretic al pr edictions. Our results 
cast doubt on lRepolust et alJ J2004 and Markova et al. (2004) 
tentative suggestion to lower supergiant mass loss rates in or- 
der to reunify the empirical WLRs with the theoretical one. 
We have indeed shown that a correction for clumping is re- 
quired even for dwarfs that are expected to show Ha in ab- 
sorption (or partially filled in by wind emission). This, in turn, 
implies that the mass loss rates should be lowered, leaving 
thus a di fference between the supergiant and dwarf WLRs. 
Moreover, Repolust et alJd2004h mention that the blue Balmer 
lines in stars with Ha emission show too much wind emis- 
sion in their cores, pointing to too high Ha mass loss rates. 
Clumps at the base of the wind might account for this be- 
haviour. We finally note that the existence of a relatively deep- 
seated clumpe d region is compat i ble with recent, theor etical 
results (e.g. Feld meier et alJll997HOwocki & Pulslll999l and 
references therein) . Although the stellar sample is limited, we 
emphasized that we found evidence of wind clu mping in all 
early O dwarfs that we have analyzed so far, see B ouret et"all 
( 2003). A unique WLR for supergiants and dwarfs might even- 
tually be empirically derived with clumped wind models once 
larger stellar samples have been analyzed. We found, however, 
that wind clumping results in revising down the mass loss rates 
by a factor larger th an the factor of about 2.3 advocated by 
lRenolustetal]ll2004l). 

6. Conclusions 

We have performed a quantitative analysis of the FUSE and 
IUE spectra of two Galactic 04-type stars, one supergiant 
(HD 190429 A) and one dwarf (HD 96715), investigating the 
role of wind clumping and its resulting effect on stellar pa- 
ramet ers and thus extending our recent st udies of SMC O 
stars jBouret et alJl2003l iHillier et alJl2003l) . Our analysis is 
based on the NLTE model atmosphere programs Tlusty and 
CMFGEN, that account for NLTE metal line blanketing and 
thus provide a detailed description of the photosphere and wind 
of O stars. We have first determined the stellar and wind param- 
eters using homogeneous wind models, but failed to reproduce 
several key lines, e.g. O V/11371, like all previous analyses of O 
stars. Similarly to the case of MC stars JCrowther et all20 02b: 
Hilli er et alJ l2003). we have to adopt a very low phosphorus 
abundance in order to match the strength of the P v/l/11118- 
1128 resonance lines with homogeneous wind models. We 
have then considered clumped wind models in order to test their 
ability to improve the fit to the observed spectra. 

The clumped wind models consistently improve the 
match to lines of different species, especially Pv/L11118- 
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1128, Ovil371 and Nrv/11718. The fit to strong Fe photo- 
spheric lines (similarly to AV 83, Hillier et al. 2003) and to 
Si IV/LU393-1402 (though still not excellent) is also improved 
in HD 190429 A. In both stars, we need to adopt a highly 
clumped wind model, f m = 0.04 (HD 190429 A) and /«, = 0.02 
(HD 96715), in order to match these lines. Based on m ea- 
sured phosphorus ionization fractions, iMassa et alJ J2003I) ar- 
gued that mass loss rates of O stars might be lower by a fac- 
tor of 5, which is thus cons isten t with f x < 0.1. I n agree- 
ment with lHillier et alJ d2003l) and lBouret et alJJ2003h . clump- 
ing needs to start deep in the wind, just above the sonic point, 
at velocities as low as v c \ « 30kms _1 , to reproduce the steep 
transition between the emission and absorption components of 
clu mp-sensitive lines. T his finding also supports the conjecture 
ofRepolustetaO§004) that their mass loss rates based on He- 
were too high for supergiants due to clumping. We showed that 
the basic physical reason behind the better fits is the increased 
recombination in clumps that lowers the wind ionization. We 
argue that our clumped wind models accurately predict the ion- 
ization structure of O-type star winds because of their ability 
to consistently match lines of different species and ionization 
stages. Our results, therefore, provide a considerably more ro- 
bust conclusion about the clumped nature of the wind of O stars 
than previous studies, and are in agreement with our earlier r e- 
sults on SMC O stars JHillier et all2003UBouret et al.l2 003). 

The main consequence of wind clumping is that mass loss 
rates need to be revised down substantially, here by a factor of 3 
(HD 190429A) to 7 (HD 96715). This factor is of the same or- 
der as the correction derived for the SMC O stars dBouret et alJ 
l2003h . We show that such a drastic correction remains in good 
agreement with Ha observations. Most earlier studies of O stars 
ignored the effect of deviations from standard winds models, 
i. e. assuming a globally stationary wind with a smooth den- 
sity/velocity stratification, to determine the properties of stel- 
lar winds. It was argued that "this standard analysis yields re- 
liable average models of the stellar wind" ( Kudrit zki & Pulsl 
2000), although it was recognized that such models are inher- 
ently incapable of reproducing the spectral signatures that sug- 
gest the existence of extensive wind structures. Although this 
argument may hold for deriving the wind velocity structure, 
we have now demonstrated that this is not the case for deter- 
mining mass loss rates. The present work, together with our 
recent work on SMC O stars, establishes that clumping likely 
is a general property of O star winds. Accounting for clumping 
will lead to a systematic and significant downward revision of 
mass loss rates. Because mass loss is a crucial aspect of mas- 
sive star evolution, a revision of mass loss rates accounting for 
the effect of clumping is urgently needed. We note in paral- 
lel that two analyses of low-luminosity O stars also revealed 
mass loss r ates much lower than those p redicted by the theoret- 
ical WLR dBouret et alJl2003t iMartins et alJl2004h . Our study 
therefore calls for a fundamental revision in our understanding 
of mass loss and of O-type star stellar winds. 

Finally, the surface abundance of HD 190429 A is in agree- 
ment with its advanced evolutionary stage, in particular show- 
ing CNO-cycle processed material (nitrogen overabundance, 
carbon and oxygen depletion) at the stellar surface. Similarly 
to SMC O dwarfs, we find also enhanced nitrogen at the sur- 



face of HD 96715, though the carbon depletion is much milder 
and oxygen is close to the solar abundance. 
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